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SECTION  I 


INTRODUCTION 


The  bulk  of  information  available  regarding  the  characteristics 
and  performance  of  multiplex  (MUX)  bus  systems  is  empirical  in 
nature.  Consequently,  designers  and  system  planners  have  had  to 
hedge  when  asked  to  estimate  voltage  waveforms,  driving  point 
impedances  and  power  distribution  for  a proposed  multiplex  bus 
system.  Either  educated  guesses  based  on  past  experience  are  made 
or  a laboratory  duplicate  of  the  candidate  system  is  built  on  which 
the  required  parameter  measurements  can  be  made.  When  relying  on 
educated  guesses  the  possibility  of  error  has  to  be  tolerated, 
especially  when  the  system  under  consideration  is  radically 
different  from  the  majority  of  existing  designs.  While  measurements 
on  a laboratory  system  generally  give  more  trustworthy  results,  they 
may  require  considerable  time  and  expense.  This  paper  describes  a 
third  alternative,  computer  simulation,  which  offers  a reasonable 
compromise  between  the  two  extremes. 

The  use  of  computer  simulation  rather  than  physical  realization 
offers  several  significant  advantages.  First,  it  greatly  simplifies 
system  analysis  when  the  bus  topology  must  be  changed  to  investigate 
various  cable  routing  schemes.  In  the  initial  design  phase  when  the 
optimum  bus  arrangement  and  routing  is  usually  uncertain,  the  use  of 
simulation  could  eliminate  the  need  for  costly  laboratory  mock-ups. 
Also,  a knowledge  of  the  system's  sensitivity  to  variations  in  cable 
or  component  parameters  is  of  considerable  importance  but  difficult 
to  obtain  in  the  laboratory.  Finally,  a need  for  worst-case 
information  in  many  cases  leads  to  the  application  of  physical  "cut 
and  try"  analysis.  If  applied  to  a complex  MUX  Bus,  this  method 
would  likely  be  excessively  expensive  and  time  consuming. 

The  computer  simulation  described  in  this  paper  may  be  applied 
to  any  MUX  Bus  which  meets  the  architectural  constraints  of  MIL-STD- 
1553  (1)  (see  Section  III).  Also,  because  of  the  computer  program's 
reliance  on  steady-state  analysis,  systems  which  exhibit  significant 
transient  response  cannot  be  completely  simulated.  However, 
laboratory  investigation  has  shown  that  the  transient  response  is 
not  significant  when  the  bus  architecture  of  MIL-STD-1553  is  used. 
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SECTION  II 


CONCLUSIONS 


As  a result  of  this  investigation,  several  observations  of 
potential  importance  to  MUX  Bus  system  planners  have  been  made. 

1)  A digital  computer  steady-state  simulation  of  MIL-STD-1553 
compatible  multiplex  buses  has  been  developed  and  tested 
against  laboratory  measurements,  yielding  acceptable 
levels  of  agreement.  This  simulation  provides  driving 
point  impedances,  power  distributions  and  voltage 
waveforms  at  various  locations  throughout  the  bus.  The 
ability  to  change  both  system  architecture  and  component 
parameters  in  software  makes  this  simulation  a valuable 
tool  for  analysis  of  proposed  MUX  Bus  configurations. 

2)  Laboratory  observations  and  simulation  results  suggest 
that  an  engineering  trade-off  must  be  made  between  power 
transferred  from  transmitter  to  receiver  along  the  bus  and 
the  level  of  signal  distortion  in  the  received  waveform. 
Power  transfer  is  greatest  when  the  stub  driving  point 
impedances  are  low.  However,  low  stub  impedances  increase 
the  loading  on  tne  bus  thereby  causing  an  increase  in 
waveform  distortion.  Computer  simulation  may  help  to  make 
this  trade-off  a more  rational  one. 

3)  It  has  been  found  that  under  certain  situations  normal 
manufacturing  variation  of  component  parameters  may 
drastically  alter  bus  performance.  The  only  tractable 
method  available  for  determining  these  sensitivities 
appears  to  be  computer  simulation. 
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SECTION  III 


MIL-STD-1553  (USAF) 


In  an  effort  to  define  operational  characteristics  and  avionics 
interfaces  for  multiplex  buses  used  in  intra-aircraft  signaling;  the 
Air  Force  Aeronautical  Systems  Division  (AFASD)  has  published  a 
document  entitled  "Military  Standard  for  Aircraft  Multiplex  Data 
Bus,"  (MIL-STD-1553  (USAF)).  The  following  is  a listing  of 
characteristics  defined  in  the  standard  which  are  relevant  to  tnis 
discussion. 


Transmission  Medium: 

Bus  Discipline: 

Transmission  Bit  Rate: 
Modulation : 

Tirri  net : 

Nurrber  of  Terminals: 

1UX  Bus  Length: 

Length  of  Terminal  Stubs: 
Short  Circuit  Protection: 

Stub  Driving  Point  Impedance 
(Receive  Mode): 


Twisted,  shielded  pair, 

Zo  = 71  ohms 

Command /Response , half 
duplex;  bus  control  exer- 
cised by  a Bus  Control  Unit 

1 Megabit  per  second 

PCM,  using  Manchester  II  code 

Asynchronous,  carried  by 
Manchester  Code 

33  maximum  including  Bus 
Control  Unit 

300  ft.  max. 

20  ft.  max. 

54  ohms  resistor  in  each 
conductor  of  the  stub  pair. 

> 6800  ohms  dc  to  100  kHz 
>2000  ohms  100  kHz  to  1 MHz 


Figure  1 depicts  the  bus  architecture,  defined  by  MIL-STD-1553, 
which  was  used  as  a basis  for  the  simulation.  For  the  cases 
examined  in  this  paper,  the  following  values  were  chosen  for  the 
variables  of  Figure  1: 

Bus  Length:  L(bus)  = 250' 

Number  of  Stubs:  N = 20 
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Length  of  each  Stub: 
Transformer  Turns  Ratios: 


Isolation  Resistances : 


Receiver  Input  Impedance: 


L(s  tub ) = 20' 

N1:N2  1 1:1  Pulse 

Engineering , 
Inc. 

N 3 : N 4 = 1:1  Type  5163 

R1  = 54  ohms 
R2  = 0 ohms 
6600  ohms 


To  assure  that  reasonable  values  were  chosen  for  the  cable 
characteristics  employed  in  the  simulation,  Haveg  Industries  type 
LE-572-0003/0002  cable  was  used  as  a guide.  This  cable,  a 
lightweight  version  of  RG-108A,  has  been  recommended  for  use  in 
multiplex  bus  applications  on  the  B-1  aircraft  by  North  American 
Rockwell.  Laboratory  measurements  on  samples  of  this  cable 
produced  the  following  parameter  values: 

C = 21.8  pf  per  foot, 

L = . 1 09  y h per  f oot , 

R = .0288  ohms  per  foot, 

Zo  = 7 1 ohms 
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SECTION  IV 


CELL  APPROACH 


definition 


In  the  MUX  Bus  system  of  MIL-STD- 1 553 , each  Multiplex  Bus 
remote  Terminal  Unit  (RTU,  see  Figure  1)  is  capable  of  transmitting 
information  onto  the  bus.  Of  course,  only  one  may  transmit  at  a 
t ime . If  the  computer  simulation  is  to  be  sufficiently  flexible  to 
allow  this  type  of  operation,  some  form  of  modularization  is 
required.  The  "T"  section  of  Figure  1 will  be  defined  as  a basic 
module,  or  "cell” , of  the  rlUX  Bus.  Each  cell  may  contain  one  RTU, 
two  transformers,  four  fault  isolation  resistors  and  three  sections 
( f shielded,  twisted  pair  cable.  In  addition,  the  first  (leftmost) 
and  last  (rightmost)  cells  each  contain  a bus  termination  resistor 
whose  value  is  equal  to  the  cable's  characteristic  impedance,  Zo. 
j ne  division  point  between  adjacent  cells  is  chosen  such  that  only  N 
cells,  where  N is  the  number  of  stubs,  are  required  to  represent  a 
riven  MUX  Bus. 

Circuit  Analysis  of  a Cell 

Transmitting  Cell 

Consider  the  cell  shown  in  Figure  2 whose  RTU  is  transmitting 
onto  the  bus.  Complex  impedances  Z(1)  and  Z(2),  the  driving  point 
impedances  seen  looking  into  the  adjacent  cells,  will  be  assumed 
known  for  the  moment.  _The  firs_t^step  taken  by  the  simulation  is  to 
determine  values  for  Z(4)  and  Z(5)  by  use  of  the  following 
relationships: (3) 

Z~(5)  = Z(2)  Cosh(yL2)  + Zo  SinhCrL2) 

Zo  Zo  Cosh(y“-L2)+  Z~(l)  SinhCT'L2) 

Z*(4)  ZO)  CoshCy'*Ll)+  Zo  SintiC^-Ll) 

IZ (2) 

Zo  Zo  CoshCy*Ll)+  Z (1)  SinhOrLl) 


where:  Y = propagation  constant  of  the  cable 

Zo  = characteristic  impedance  of  the  cable 

Impedances  Z(6),  Z(6),  Z(3),  Z(9),  and  Z(7)  are  obtained  in  a 
similar  manner  using  the  above  relationship  combined  with  steady 
state  circuit  analysis  and  the  transformer  lumped  component 
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equivalent  circuit.  This  procedure  must  be  repeated  at  each 
frequency  which  corresponds  to  an  harmonic  component  of  the  RTU 
output,  voltage  V(7)  of  Figure  2.  The  waveform  in  Figure  3 is  a 
normalized  representation  of  this  unit's  output  voltage.  The 
Fourier  series  for  this  voltage  is: 


V (7) 


E 

ns=  1 


C cos(nwt) 
n 


(3) 


where: 


Cn 


„ . ntl 

Sin  tt 

[sin  tt  n(to  ± U)1 
T 

ntl 

r x J 

n(to  + tl) 

TT  

L T 

n = harmonic  number  = 1,  3,  5, 

Therefore,  as  many  odd  harmonies  of  the  fundamental  frequency  as 
accuracy  requirements  dictate  will  have  to  be  considered  in  the 
impedance  calculations.  [NOTE:  For  rise/fall  times  consistant  with 

MI L— STD— 1 553  and  a 1 MHz  fundamental  frequency,  approximately  six  or 
more  odd  harmonies  should  be  considered.] 

After  the  impedances  Z(1)  through  Z(9)  are  known  at  each 
applicable  frequency,  the  next  step  is  to  apply^each  harmonic 
component  of  the  JTU  outpqt  voltage  ^wavef orm  V(7)  to  the  cell. 
Calculation  of  V(9)  and  V ( 3 ) from  V ( 7 ) is  fairly  straightforward  and 
will  not  _be  explained  in  detail.  Voltage  V ( 8 ) is  found  by  using 
voltage  V(3)  and  the  following  eauation  (3): 


V(8)  = V(3)  [CoshOy • L3) - ( Z0/ZO ))  SinhCg-L3)]  (4) 


Using  this  relationship,  the  transformer  model,  and  AC  circuit 
analysis,  voltages  V(6),  V(4),  V(5),  V ( 1 ) and  V(2)  may  be  obtained. 
As  with  the  impedance  calculations,  this  procedure  must  be  repeated 
for  each  of  the  relevant  harmonic  frequencies. 

The  final  step  in  the  analysis  of  the  transmitting  cell  is  to 
determine  the  composite  voltage  waveforms,  harmonic  power 
distributions,  and  composite  power  distributions  at  each  of  the  nine 
points  shown  in  Figure  2.  The  composite  voltage  waveform  at  point  Q 
is  found  by: 

v(t)  = V'  V^(Q)  Cos(nwt)  (5) 

Q (n=lTf,5...) 

[composite  voltage] 
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where: 


n = harmonic  number 


NN  = highest  odd  harmonic  considered 
Q = position  number  (1,2,..., 9) 

Vn(Q)  = voltage  phasor  at  point  Q for  frequency  n 

In  the  computer  program  this  summation  is  done  at  discrete  time 
increments  from  t=0  to  t=T,  where  T is  the  period  of  the  fundamental 
frequency. 

The  impedances  previously  calculated  are  actually  driving  point 
impedances  at  each  of  the  measurement  points  indicated  in  Figure  2. 
Using  this  fact,  along  with  the  corresponding  harmonic  voltages, 
leads  to  the  following  expression  for  average  harmonic  power  at 
point  Q: 


MQ) 


lV(Q)l2 

2 | Zn (Q)  | 


cos  (0) 


average 

harmonic 

power 


(6) 


where:  0 = angle  of  impedance  Z^Q) 

To  find  the  composite  average  power  at  point  Q requires  a 
summation  of  the  average  harmonic  powers,  i.e.: 


P 


Q 


2 


pn(Q) 


composite 

average 

power 


(7) 


The  only  matter  not  covered  at  Jbhis  poiqt  is  where  values  for 
the  previously  assumed  impedances  Z(1)  and  Z(2)  come  from.  This 
question,  addressed  in  the  next  two  sections,  is  the  basis  for  the 
flexibility  inherent  in  the  cell  approach. 

Receiving  Cell 

Within  the  framework  of  MIL-STD- 1 553 , only  one  RTU  can  transmit 
information  onto  the  bus  at  any  given  time.  Consequently,  (N-1) 
cells  (where  N is  the  number  of  RTU's)  are  considered  to  be 
receiving  cells.  For  the  purpose  of  this  discussion  these  (N-1) 
cells  will  be  defined  as  either  Left  Cells  (LCELL)  or  Right  Cells 
(RCELL),  depending  on  their  location  relative  to  the  transmitting 
cell. 
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Right  Cells.  Figure  4 depicts  impedance  and  voltage 
conventions  within  an  RCELL.  These  conventions  are  applied  to  all 
RCELL's  except  the  last , or  rightmost  cell.  The  difference  in 
RCELL(LAST)  is  that  f{2)  is  actually  equal  to  RTERM,  a known 
constant,  as  can  be  seen  from  Figure  1.  The  significance  of  this 
point  will  be  discussed  in  the  section  on  chaining  of  cells. 

As  in  the  transmitting  cell  case,  the  first  step  in  the 
simulation  of  an  RCELL  is  to  determine  the  driving_j)oint  impedances 
throughout  the  cell.  Assuming  that  the  value  of  Z ( 2 ) is  known: 


Z(5)  = Z~(2)  Cosh(T*L2)+  Zo  Sinh(yL2) 
Zo  Zo  Cosh(yL2)  + ZU)  Sinh(yL2) 


(8) 


Mote  that  this  equation  has  the  same  form  as  Equation  1. 

The  receiving  input  impedance  of  the  RTU,  zT7),  is  known. 
Therefore,  impedances  Z(9)  and  Z ( 3 ) can  be  readily  calculated.  To 
find  Z(8),  zt*3)  is  used  in  the  following  equation: 


Z (8)  _ Z (3 ) CoshCyL3)+  Zo  SinhCyL3) 
Zo  • Zo  Cosh(y*L3)+  Z?3)  Sinh(Y*L3) 


Solving  for  Z(6)  using  the  transformer  model  and  Z(8)  allows 
the  following  relation  for  Z(4): 


Z (4)  = [z (6)  + 2-Rl] | | Z (5)  (10) 

With  the  value  of  Z ( 4 ) determined,  Z ( 1 ) can  be  found  by: 

Z(l)  = Z~f4)  Cosh(y-Ll)+  Zo  SinhCy-Ll) 

Zo  Zo  Cosh(y*Ll)+  Z(4)  SinhCyLl) 

An  important  point  to  be  made  here  is  that  ZO)  of  this  cell  is 
identically  equal  to  Z(2)  of  the  preceding  cell  on  the  left. 
Conversely,  Z(2)  is  equal  to  zfl ) of  the  following  cell  on  the 
right.  This  transference  of  impedance  values  across  cell  boundaries 
provides  a "daisy-chaining"  effect,  supplying  values  for  Z72)  of  the 
transmitter  cell,  which  was  assumed  known  in  the  previous 
discussion. 

The  next  step  in  the  analysis  is  to  determine  the  voltages 
within  the  cell.  Here  the  input  voltage  is  considered  to  be  V(1), 
as  opposed  to  vT7)  in  the  transmitting  cell  case.  Values  for  this 
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voltage  at  the  various  harmonics  are  supplied  by  V(2)  of  the  cell  on 
the  immediate  left.  To  find  V(4)  apply  vTT)  to: 

V(4)  = vfl)  [CoshCy ’Ll) -(Zo/zTl) ) SinhCy-Ll)]  (12) 


which  has  the  same  form  as  Equation  (2).  Using  the  voltage  divider 
relationship: 


V (6)  = V (4) 


_ Z(6) 
Z(6)  + 2*Rl 


(13) 


Solving  for  V(d)  using  the  transformer  model  leads  to: 


V (3)  = V(8)  [CoshCy-L3)-  (Zo/Z?8))  Sinht?-L3)]  (14) 


Calculating  V ( 9 ) from  the  transformer  model,  determine  the  received 
harmonic  voltage  at  the  RTU  by: 


V (7)  = V (9) 


Z(7) 

Zf7)  + 2*R2 


where:  Z(7)  represents  the  RTU  receiver  input  impedance. 


(15) 


Returning  to  the  bus/stub  junction,  apply  the  already  determined 
value  for  vT^)  to: 


V(2)  = V (4)  Coshty*L2)-  (Zo/zX5))  Sinhty*L2) 

As  mentioned  earlier,  voltage  V(2)  is  the  input  voltage  V(1)  for  the 
next  RCELL  to  the  right. 

After  the  above  procedures  have  been  used  to  obtain  the 
impedances  and  voltages  for  each  harmonic,  Equations  5,  6,  and  7 are 
used  to  determine  the  composite  voltage  waveforms,  average  harmonic 
powers  and  average  composite  powers  in  each  of  the  RCELLS. 

Left  Cells.  The  analysis  of  LCELLS,  cells  which  are  to  the 
left  of  the  transmitting  cell,  is  analogous  to  RCELLS.  The  voltage 
and  impedance  conventions  for  this  case  are  shown  in  Figure  5. 

Again,  the  difference  in  the  first  or  leftmost  cell  is  that  zTl ) is 
equal  to  RTERM,  a point  that  will  be  discussed  further  in  the  next 
section. 
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In  an  LCELL  the  entry  point  for  impedance  calculations  is  Z(1). 
Values for  this  impedance  at  each  harmonic  frequency  are  obtained 
from  Z(2)  in  the  previous  cell  to  the  left.  The  order  of 
calculation  and  methods  used  for  _an  LCELL  are  the  same  as  for  an 
RCELL,  i.e.,  Z ( 4 ) is  found  from  Z ( 1 ) and  the  ^impedance  translation 
equation;  impedances  Z(9),  Z(3),  Z(8),  and  Z ( 6 ) are  obtained  by 
using  the  transformer  model,  translation  equation,  and  AC  circuit 
analysis;  Zj2)  is  determined  by  using  the  parallel  combination  of 
Z ( 4 ) and  Z(5)  in  the  translation  equation.  Calculation  of  voltages 
within  an  LCELL^  is  also  analogous  to  RCELL's.  Instead  of  using  V ( 1 ) 
as  the  input  V(2)  is  used,  but  the  procedure  to  be  followed  is 
equivalent.  The  main  difference  between  the  analysis  of  an  LCELL 
and  an  RCELL  is  the  direction  of  the  drivinr  point  impedances,  and 
the  variables  which  are  used  as  inputs  to  the  cell.  This  similarity 
greatly  decreases  the  complexity  of  the  computer  simulation  program. 
Also,  it  allows  complete  flexibility  of  transmitter  location--a  fact 
that  will  be  illustrated  in  the  following  section. 

Chaining  of  Cells 

Having  considered  each  of  the  different  cell  types  separately, 
the  next  step  is  to  show  how  they  may  be  connected  to  form  a 
complete  multiplex  bus.  Consider  the  system  of  Figure  6 in  which  a 
transmitting  cell,  LCELL,  and  RCELL  have  combined  to  form  a three 
subscriber  MUX  Bus. 

The  first  step  in  the  analysis  of  this  system  is  calulation  of 
the  driving  point  impedances.  Referring  to  Cell  3 of  Figure  6,  an 
RCELL,  several  points  about  the  cell  boundaries  become  apparent.  In 
the  discussion  of  impedance  calculations  in  an  RCELL,  mention  was 
made  that  the  value  of  Z ( 2 ) would  be  supplied  by  zTl ) of  the  next 
RCELL  to  the  right.  However,  in  this  case  Zl2)  is  actually  equal  to 
ZTERM,  a known  constant.  Therefore,  impedances  zTl ) through  zT9) 
can  be  calculated  for  this  cell  without  regard  to  any  other  cell. 

If^ there  were  more  than  one  RCELL,  the  value  just  calculated  for 
Z ( 1 ) in  the  rightmost  cell  would  become  the  ZT2)  of  the  next  RCELL 
to  the  left,  with  the  process  continuing  until  all  of  the  RCELLS  had 
been  considered.  The  same  procedure  is  applied  to  the  LCELLs  except 
that  calculation  begins  at  the  leftmost  cell  in  which  zTl)  = ZTERM 
and  then  progresses  to  the  right  until  all  LCELLS  have  been 
addressed.  At  that  point  the  only  driving  point  impedances  not  yet 
calculated  would  be  those  in  the  transmitting  cell  itself.  Since 
there  is  no  requirement  that  the  number  of  LCELLS  be  equal  to  that 
of  RCELLS  the  location  of  this  remaining  transmitting  cell  is 
completely  arbitrary,  a major  advantage  of  the  cell  approach. 

Returning  to  the  transmitting  cell  of  Figure  6,  having  found 
values  for  Z ( 2 ) of  the  LCELL  on  the  immediate  left  and  Z~0 ) of  the 
RCELL  on  the  immediate  right,  calculation  of  the  driving  point 
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Figure  6 THREE  CELL  MULTIPLEX  BUS 


impedances  within  the  transmitting  cell  follows  the  discussion  in 
the  Transmitting  Cell  Section. 

Calculation  of  voltages  throughout  the  bus  proceeds  in  an 
analogous  fashion  but  with  an  opposite  direction  of  flow.  Here  the 
known  values  are  the  harmonic  voltage  amplitudes  of  the  transmitting 
RTU's  output  waveform.  Since  all  impedances  have  been  determined, 
it  is  a simple  matter  to  find  the  voltages  within  the  transmitting 
cell  using  these  amplitudes  and  the  methods ^described  previously. 

The  values  produced  for  voltages  and  V(2)  become  the  inputs  for 

the  cells  on  the  immediate  left  and  right.  The  chaining  process 
continues  until  all  cell  voltages  on  each  side  of  the  transmitting 
cell  have  been  calculated.  This  sequence  of  operations  may  be 
summarized  as  follows:  impedance  calculations  flow  from  the 

extremities  of  the  mux-bus  inward  toward  the  transmitting  RTU,  while 
voltage  calculations  flow  from  the  transmitting  RTU  outward  toward 
the  extremities  of  the  bus. 

With  all  voltages  and  impedances  determined  throughout  the  bus, 
calculation  of  the  composite  voltage  waveforms,  average  harmonic 
powers,  and  average  composite  powers  for  each  cell  proceed  as 
previously  described  in  the  Transmitting  Cell  Section. 
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SECTION  V 


MODEL  VALIDATION 


Validation  of  the  MUX  Bus  simulation  was  accomplished  in  two 
stages.  First,  hand  calculation  and  Smith  chart  analysis  were  used 
to  check  the  driving  point  impedance  and  power  distribution 
algorithms.  Next,  the  voltage  waveforms  produced  by  the  simulation 
were  compared  with  oscillographs  taken  on  a laboratory  mock-up  of  a 
representative  MUX  Bus  (2).  The  bus  architecture  used  for  this 
comparison,  shown  in  Figure  7,  had  the  characteristics  discussed  in 
Section  III  with  one  exception.  Since  it  was  possible  to  obtain  only 
a limited  supply  of  pulse  transformers,  the  staggered  transformer 
placement  arrangement  shown  was  adopted. 

CASE  I.  The  first  comparison  between  laboratory  waveforms  and 
computer  simulated  plots  is  shown  in  Figure  3.  The  transmitter,  a 
square  wave  pulse  generator  with  output  impedance  equal  to  Zo,  was 
connected  to  the  end  of  Stub  6.  The  point  at  which  voltage 
waveforms  were  to  be  compared  was  the  junction  between  Stub  4 and 
the  bus  (labeled  point  "A"  on  Figure  7).  The  laboratory 
oscillograph  of  this  waveform  when  isolation  resistors  HI  (Figure  1) 
were  included  in  all  cells  is  shown  on  the  left  of  Figure  8.  The 
dashed  line  on  the  plot  directly  above  this  picture  is  the  computer 
simulation  of  the  same  voltage  measured  over  one  time  period  and 
normalised  to  a +1  volt  transmitter  signal  excursion.  The  solid 
line  in  tnis  plot  is  the  simulation  plotted  waveform  for  the 
reference  point  defined  by  MIL-STD-1553  (sec  point  "R",  Figure  7). 
For  the  sake  of  clarity,  the  oscillograph  of  this  reference  waveform 
is  not  snown. 

The  right  half  of  Figure  8 shows  the  comparison  between 
laboratory  and  simulated  waveforms  for  the  same  measurement  point 
(point  "A”,  Figure  7)  when  the  isolation  resistors  R1  are  shorted 
out.  Again,  the  dashed  line  on  the  computer  generated  plot 
represents  the  simulations  approximation  of  the  waveform  at  this 
point. 

When  each  set  of  comparisons  is  studied,  two  facts  become 
apparent.  First,  the  computer  simulation  plotted  waveforms  agree 
very  closely  with  their  associated  laboratory  oscillographs.  This 
implies  that  the  simulation  is  indeed  valid.  The  second  observation 
involves  the  utility  of  the  isolation  resistors  (R1).  As  shown  on 
the  left  of  Figure  8,  the  observed  and  simulated  waveforms  are 
relatively  free  from  distortion  when  R1  is  included  in  each  stub. 
However,  as  shown  on  the  right,  when  this  resistance  is  removed, 
there  is  a significant  degree  of  "ringing"  apparent.  This 
distortion  is  a measure  of  the  mutual  interaction  between  stubs. 
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Figure  7 TWENTY  CELL  MULTIPLEX  BUS 
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Since  the  degree  of  interaction  is  an  inverse  function  of  the 
driving  point  impedance  of  the  stubs,  the  "smoothing"  caused  by  the 
inclusion  of  resistances  R1  is  understandable. 

CASE  II.  The  second  comparison  to  be  discussed  refers  to  the 
same  bus  architecture  as  CASE  I,  but  considers  the  waveform  at  point 
"B"  of  Figure  7*  The  laboratory  oscillographs  and  computer 
simulated  waveforms  with  and  without  isolation  resistances  (R1)  are 
shown  in  Figure  9.  As  in  the  previous  case,  good  agreement  between 
the  oscillographs  and  the  corresponding  dashed  waveforms  on  the 
computer  plots  empirically  justifies  the  simulation  algorithm.  In 
addition,  the  set  of  waveforms  for  the  R1  = 0 situation  graphically 
displays  the  type  of  problem  that  can  occur  if  a high  degree  of 
interaction  is  allowed  between  stubs.  If  a zero-crossing  detector 
were  used  on  this  waveform,  an  erroneous  output  would  probably  occur 
if  la;  pass  filtering  were  not  used.  However,  since  the  simulation 
can  predict  tnis  potential  problem  the  designer  could  employ 
filtering,  isolation  resistors,  or  a different  type  of  detector  to 
avoid  this  situation  in  an  operational  system. 

CASE  III.  The  final  comparison  to  be  made  is  shown  in  Figure 
10.  As  in  the  previous  cases,  the  bus  architecture  and  transmitter 
location  are  unchanged.  Here  the  point  of  observation  has  been 
moved  to  "C"  in  Figure  7.  Good  agreement  is  again  shown  but,  more 
importantly,  the  simulation  exhibits  the  same  bandwidth  limita- 
tion on  a long  bus  that  was  found  during  the  laboratory  investiga- 
tion . 


By  considering  the  cases  discussed  above,  along  with  many 
others  not  shown,  a reasonable  degree  of  confidence  in  the  accuracy 
of  the  simulation  was  obtained.  Hence,  the  next  step  in  the 
investigation  was  to  exercise  the  simulation  beyond  the  scope  of  the 
laboratory  effort  in  an  attempt  to  expose  possible  problems  or 
improvements  in  MIL-STD- 1 553-  One  of  the  early  results  of  this 
study  is  shown  in  Figure  11. 

The  question  addressed  here  was  the  effect,  if  any,  on  voltage 
waveforms  caused  by  variation  in  the  inductance  of  the  cable  used 
for  the  bus.  A laooratory  analysis  of  this  question  would  be 
impractical,  to  say  the  least.  However,  as  shown  by  the  plotted 
waveforms  in  Figure  11,  a possible  problem  area  would  be  missed  if 
this  study  were  not  done.  In  all  four  plots  the  dashed  waveform  is 
the  simulated  voltage  at  point  "B"  of  Figure  7.  Case  II  previously 
discussed  this  voltage  waveform  for  a nominal  cable  inductance  of 
0.109  /ih  per  foot.  Plots  "A"  and  "B"  on  the  left  of  Figure  11  show 
the  waveforms  produced  by  varying  the  cable  inductance  +.10^  about 
its  nominal  value  when  isolation  resistances  R1  were  installed  in 
the  bus  of  Figure  7.  Comparing  these  two  plots  it  is  apparent  that 
little  sensitivity  to  an  inductance  variation  of  this  magnitude 
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would  be  expected.  However,  consider  plots  "C"  and  "D"  on  the  right 
of  Figure  11.  Here  a considerable  change  in  waveform  is  displayed 
when  resistors  R1  are  removed  from  the  stubs  and  L is  varied  over 
the  same  ±*\0%  range.  Since  cable  manufacturing  tolerances  are 
within  the  range  considered,  it  is  entirely  possible  that  a designer 
expecting  waveform  "D"  might  instead  find  waveform  "C"  when  the  bus 
was  actually  constructed,  perhaps  causing  incompatibility  with  his 
detection  circuitry.  This  inherent  ability  of  computer  simulation 
to  permit  variation  of  any  system  parameter  by  a simple  change  in 
software  provides  a powerful  tool  for  analysis  of  multiplex  bus 
systems. 
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APPENDIX  I 


SIMULATION  SOFTWARE 


In  order  to  provide  a better  understanding  of  the  simulation 
and  to  facilitate  its  use  by  interested  parties,  the  following  four 
topics  will  be  addressed: 


Part 

1 - 

The  computer  system  on  which  the  simulation  was 
exercised,  including  required  capacities  and 
peripheral  equipment. 

Part 

2 - 

The  basic  simulation  routine,  called  TPM0D2,  with  a 
discussion  of  required  inputs  and  available  outputs. 

Part 

3 - 

The  plotting  routine,  PLM0D2,  which  produced  the 
voltage  waveform  plots  shown  in  Section  IV. 

Part 

4 - 

The  operational  steps  required  to  exercise  the 
simulation  on  a computer  installation  similar  to  the 
author's. 

Part 

1 : 

Computer  System  Requirements  - The  simulation  actually 

consists  of  two  separate  programs,  TPM0D2  and  PLM0D2.  Both  of  these 
programs  were  written  in  FORTRAN  IV  and  run  on  an  IBM  370/155 
computer.  To  exercise  the  programs  without  modification,  the 
reauirements  of  Table  A1,  Figure  1.1  should  be  met.  Table  A2  lists 
the  core  and  CPU  demands  found  during  a representative  run  made  for 
a 32  stub  MUX  Bus. 

Part  2:  Description  of  TPM0D2  - In  order  to  reduce  computer 

core  requirements  and  to  increase  flexability,  the  simulation  was 
broken  into  two  parts.  The  first  part,  called  !,TPM0D2",  uses  the 
architectural  and  component  parameter  inputs  to  produce  a magnetic 
tape  on  which  the  complex  voltage  "phasors"  for  each  harmonic  at 
each  node  of  Figure  1 are  written.  The  operation  of  this  program 
can  best  be  explained  by  reference  to  the  flow  chart  of  Figure  1.2 
and  the  listing  of  the  program  immediately  following. 

TPM0D2  has  essentially  a three  tiered  structure.  The  first 
level,  or  main  program,  handles  input  and  output,  generation  of  the 
input  voltage  harmonic  components,  power  calculations,  and  the 
sequencing  of  subroutine  calls.  Since  this  analysis  is  multi- 
variate, involving  both  positions  along  the  bus  and  frequency 
(harmonic),  there  are  several  nested  loops  in  the  sequencing  portion 
of  the  routine.  The  harmonic  number,  or  frequency,  forms  the  outer 
loop  inside  of  which  the  cell  number,  or  position  along  the  bus,  is 
varied.  As  these  parameters  take  on  their  allowed  doublets,  calls 
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COMPUTER  SYSTEM  REQUIREMENTS 

TO  EXERCISE  TPMOD2  WITHOUT  MODIFICATION 

TO  EXERCISE  PLM002  WITHOUT  MODIFICATION 

1)  FORTRAN  12  COMPILER 

2)  COMPLEX  NOTATION  a ARITHMETIC 

3)  9-TRACK  TAPE  DRIVE 

1)  FORTRAN  12  COMPILER 

2)  COMPLEX  NOTATION  8 ARITHMETIC 

3)  9-TRACK  TAPE  DRIVE 

4)  7-TRACK  TAPE  DRIVE 

5)  CALCOMP  PLOTTER  AND 

SOFTWARE  PACKAGE 

TABLE  A I 


REPRESENTATIVE  CORE  REQUIREMENTS 
AND  CPU  TIMES  FOR  A 32  STUB  MUX-BUS  SIMULATION 

STEP 

TPMOD2 

PLMOD2 

( 1 3 PLOTS  PRODUCED  ) 

COMPILE 

CPU  CORE 

I5S  I20K 

CPU  CORE 

6S  I04K 

LINK/EDIT 

2 S 1 22  K 

2S  I22K 

RUN 

7S  246K 

IIS  242  K 

TABLE  A2 


Figure  1. 1 COMPUTER  SYSTEM  REQUIREMENTS 
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Figure  1.2  FLOW  CHART  OF  TPM0D2 
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VFORM 


(VOLTAGE 
TRANSLATION 
THRU  THE 
TRANSFORMER) 


Figure  1.2  FLOW  CHART  OF  TPMOD2  ( CONT.) 
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are  made  to  the  routines  in  the  second  tier:  LCELLZ,  RCELLZ,  LVOLT, 

and  RVOLT.  Each  of  these  subroutines  controls  the  calculation 
sequence  within  a given  cell,  making  calls  to  the  four  service 
routines,  IMPDNS,  EMF,  ZFORM,  and  VFGRM  when  appropriate.  IMPDNS 
and  EMF  exercise  the  transmission  line  translational  relationships 
for  impedance  and  voltage,  respectively.  VFORM  and  ZFORM  transform 
voltages  and  impedances  through  each  of  the  two  transformers  that 
may  be  contained  in  a stub. 

The  complex  voltages  and  impedances  calculated  by  the  above 
procedure  for  each  relevant  harmonic  of  the  input  waveform  are 
stored  in  memory.  When  all  values  have  been  found,  the  voltage 
phasors  are  written  onto  a 9-track  magnetic  tape  for  later  use  by 
the  plotting  routine,  PLM0D2.  The  power  transfer  ratios  from  the 
transmitter  to  the  receivers  are  then  calculated  and  printed  along 
with  the  complex  impedances  previously  determined. 

Input  Data  Format  - The  input  data  required  by  the  program, 
including  format  and  card  column  designations,  is  shown  in  Figure 
1.3.  The  column  marked  "Data  Card  Number"  refers  to  the  order  in 
which  the  data  cards  appear  in  the  source  deck. 

Output  Formats  - In  addition  to  the  complex  voltages  written  on 
magnetic  tape,  TPM0D2  produces  hard  copy  listings  of  the  following: 

a)  The  input  data,  including  the  calculated  value  of  the 
cable  characteristic  impedance  and  the  number  of  odd 
harmonics  considered.  (See  Figure  1.4.) 

b)  The  driving  point  impedance  seen  by  the  transmitter  as  a 
function  of  frequency  (Figure  1.5). 

c)  The  driving  point  impedance  of  a detached  stub  as  seen  by 
the  main  bus  (Figure  1.6). 

d)  The  driving  point  impedance  looking  into  each  junction  of 
a stub  and  the  main  bus.  For  cellj^to  the  left  of  the 
transmitter,  this  corresponds  to ^(5)  of  Figure  5,  while 
for  cells  on  the  right  this  is  Z(4)  of  Figure  4^  In  the 
transmitter  cell  this  impedance  corresponds  to  Z(4)  in 
parallel  with  Z(5).  Also  included  is  the  power  loss  from 
the  transmitter  to  the  junction,  both  as  a function  of 
frequency  and  total  (Figure  1.7). 

e)  The  power  loss  from  the  transmitter  to  each  receiving  RTU, 
combined  and  as  a function  of  frequency.  This  data,  shown 
in  Figure  1.8,  is  an  essential  input  to  the  determination 
of  required  transmitter  power  levels. 


IA- 45,004 


DATA 

CARD 

NUMBER 

VARIABLE 

NAME 

DESCRIPTION 

CARD 

COLUMNS 

FORMAT 

COMMENTS 

1 

Nl(  1 ) 

TRANSFORMER  Tl 
PRIMARY  TURNS 

1 - 5 

F 5. 2 

N,  OF  FIG  2 

N2(l  ) 

Tl  SECONDARY 
TURNS 

6-10 

F 5 2 

N,  OF  FIG  2 

tri ( 1 ) 

Tl  PRIMARY 
RESISTANCE  (H) 

11-15 

F 5.2 

RP  OF  FIG.  2 

TR  2 ( 1 ) 

Tl  SECONDARY 

resistance 

16-20 

F 5 .2 

Rs  OF  FIG  2 

TLI  l l ) 

Tl  PRIMARY 

leakage 

inductance 

21-30 

E 10. 4 

LP  OF  FIG.  2 

Tl  2 ( l ) 

Tl  SECONDARY 

LEAKAGE 

INDUCTANCE 

31  -40 

ElO  4 

Ls  OF  FIG  2 

TLM  ( 1 ) 

Tl  MAGNETIZING 
INDUCTANCE 

41-50 

E 10.4 

ln  of  fig.  2 

2 

N I (2  ) 

T2  PRIMARY 
TURNS 

l -5 

F 5 .2 

N 2 (2) 

T2  SECONDARY 
TURNS 

6 - 10 

F 5 . 2 

TRI (2) 

T2  PRIMARY 
RESISTANCE 

H-15 

F 5 . 2 

TR2 (2) 

T2  SECONDARY 

resistance 

16  -20 

F 5. 2 

v SAME  AS 
ABOVE 

TLI (2) 

T 2 PRIMARY 

leakage  inductance 

21-30 

E i 0. 4 

TL2 ( 2 ) 

T 2 SECONDARY 

leakage  inductance 

31-40 

ElO  4 

TLM ( 2 ) 

T2  MAGNETIZING 
INDUCTANCE 

41-50 

E 10  4 

.J. 

3 

R 

CABLE  SERIES 
RESISTANCE  /FOOT 

l - 12 

EI2.6 

V 

J 

UNIT  OF  LENGTH 
IS  ARBITRARY  AS  LONG 
AS  ALL  ARE  EQUAL  , 
CHARACTER  1 ST  1 C 
^ IMPEDANCE  Zo 
IS  CALCULATED 
INTERNALLY  BY  THE 
SI  MUL  ATION 

IND 

CABLE  SERIES 

inductance/foot 

13  -24 

E 12  6 

G 

CABLE  SHUNT 
CONDUCTANCE/  FOOT 

25-36 

E I 2 . 6 

C 

CABLE  SHUNT 
CAPACITANCE  / FOOT 

37-  4B 

E 12  6 

T 

INPUT  PULSE 
WIDTH  (SECONDS) 

49  - 60 

E 12  6 

maxhar 

MAXIMUM  NUMBER 
OF  HARMONICS 
ALLOWABLE 

61-63 

I 3 

SEE  NOTE  AT  " PERCN  T *' 
ON  CARD  0 4 

Figure  1.3  TPM0D2  INPUT  DATA  CARD  FORMAT 
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I A- 45 ,085 


DATA 

CARO 

NUMBER 

VARIABLE 

NAME 

DESCRIPTION 

CARD 

COLUMNS 

FORMAT 

COMMENTS 

4 

TRISE 

INPUT  PULSE  RISE 

time 

i-»2 

E l 2. 6 

FOR  FOURIER  SERIES 
PRESENTLY  USED 
Y TRISE  = TFALL  ; 

MEASURED  FROM  0 - 1 00% 

TFALL 

INPUT  PULSE  FALL 
TIME 

13  -24 

El  2.6 

PERCNT 

PERCENT  ENERGY 
CUTOFF  FOR  INPUT 
HARMONIC  GENERATOR 

25-36 

E I 2 6 

harmonics  of  input  waveform 

ARE  GENERATED  UNTIL  ENERGY 
FALLS  BELOW  THIS  PERCENT  OF 
THE  FUNDAMENTALS  ENERGY 

2 GEN 
1 

INPUT  WAVEFORM 
GENERATORS 

37  —4  B 

EI2.6 

R«  [ZGEN  ] 

1 

INTERNAL 

IMPEDANCE 

49-60 

E 12  6 

Im  [ZGEN] 

5 

LAST 

NUMBER  OF  LAST 
CELL  ON  RIGHT 

»-2 

12 

EOUIVALENT  TO  THE 
TOTAL  NUMBER  OF  STUBS 

XMTR 

NUMBER  OF  CELL 
CONTAINING  TRANSMITTER 

3-4 

I 2 

CELLS  NUMBERED 
FROM  LEFT  TO  RIGHT 

RESI 

LEFT  BUS  TERMINATION’S 
RESISTANCE 

5 - 14 

EI0.4 

EOUIVALENT  TO  RTERM 

LI 

LEFT  BUS  TERMINATION’S 
INDUCTANCE 

15-24 

E 10.4 

GENERALLY  * 0 

rlast 

RIGHT  BUS  TERMINATION'S 
RESISTANCE 

25-  34 

E 10.4 

► SAME  AS  ABOVE 

LL  AST 

RIGHT  BUS  TERMINATION’S 
INDUCTANCE 

35  -44 

E 10.4 

J 

6-  37* 

RI(N) 

TOP  ISOLATION 
RESISTANCES 

1 - 10 

F 1 0.0 

TOTAL  TOP  ISOLATION 
RESISTANCE  = 2*  R 1 ( N ) 

R2(N) 

BOTTOM  ISOLATION 
RESISTANCES 

H-20 

FI  0 0 

RSTUB (N) 

STUB  TERMINATION'S 
RESISTIVE  COMPONENT 

rv> 

1 

O* 

o 

F 1 0 0 

DETERMINED  BY  THE 
► RECEIVING  RTU’S 
INPUT  IMPEDANCE 

LSTUB (N ) 

STUB  TERMINATION'S 
INDUCTIVE  COMPONENT 

31-40 

EI0.4 

L ( 1 , N ) 

CELL  CABLE 
length 

41-45 

15 

-S 

L(  2 , N ) 

« 

* 

<T» 

1 

LA 

O 

15 

► SEE  FIG  2 

L(3,N) 

" 

51  - 55 

15 

TXFMR (N) 

INDICATES  IF  TOP 
TRANSFORMER  IS  PRESENT 

56  -60 

15 

s 

s 1, TRANSFORMER  PRESENT 

BXFMR(N) 

INDICATES  IF  BOTTOM 
TRANSFORMER  IS  PRESENT 

61  - 65 

15 

= 0, TRANSFORMER  ABSENT 

•ONE  CARO  CONTAINING  THESE  PARAMETERS  IS  REOUIRED  FOR  EACH  OF 
THE  N CELLS  IN  THE  SYSTEM  BEING  SIMULATED 


Figure  1.3  TPM0D2  INPUT  DATA  CARD  FORMAT  (CONT.) 
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FIGURE  1.4  (Con't) 


DRIVING  POINT  IMPEDANCE  AT  XM TR/S  TUB  JUNCTION  .VS.  FRFQUENCY 
XMTR  IS  IN  CELL  14 
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**********  SERIAL  NO.  09.53. SI  ********** 


FIGURE  1.5 
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DKIVIN 


PUINT  IMPEDANCE  OF  STUB  1 .VS.  FREQUENCE 


HO  RT  Z 

n 

IHMS 

MAGNITUDE 

ANGLE 

1 OOE 

07 

934 .45 

♦ J( 

-2050.66) 

2253.53 

-65.50 

300E 

07 

254.28 

♦Jl 

-423.72 ) 

494.16 

-59.03 

500E 

07 

1 9 9. 3 p 

120.31) 

232.87 

31.11 

7000 

07 

186.10 

♦J( 

502.02 ) 

535.40 

69.66 

900E 

07 

164.06 

347.32) 

867.08 

77.74 

1 10E 

03 

190.61 

♦Jt 

1213.58  ) 

1228.46 

8 1.07 

1 30c 

OR 

218.98 

♦ J( 

1695.03 ) 

1709.1  1 

<4- 

>0 

fy 

00 

150F 

0° 

434.75 

♦ J( 

2798.05) 

2831.62 

81.17 

170E 

03 

1735.93 

-2085.90) 

2713.76 

-50.23 
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FIGURE  1.6 


55 


NUDE  DRIVING  POINT  IMPEDANCE  -VS-  FREQUENCY  FOR  CELL  1 

HCRTZ  OHMS  MAGNITUDE  ANGLE  POWER  RATIO ( DB ) 
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FIGURE 


The  information  supplied  by  these  tables,  when  combined  with 
the  waveform  plots  produced  by  PLM0D2,  provided  an  essentially 
complete  picture  of  the  performance  that  can  be  expected  from  a 
proposed  multiplexed  bus. 

Part  Description  of  PLM0D2  - The  purpose  of  PLM0D2  is  to 
reconstruct  and  plot  the  voltage  waveforms  at  the  points  considered 
by  TPM0D2.  The  operation  of  the  program  is  shown  in  the  flow  chart 
of  Figure  1.9  and  the  listing  which  immediately  follows. 

Input  Data  - Input  to  the  plotting  routine  comes  from  two 
sources.  The  9-track  magnetic  tape  produced  by  TPM0D2  provides  the 
harmonic  voltage  amplitudes  and  phases  required  to  reconstruct  the 
waveforms.  Selection  of  the  waveforms  to  be  plotted,  the  number  of 
time  increments  per  pulse  period,  and  designation  of  the  reference 
voltage  in  the  transmitter  cell  are  all  entered  as  data  cards  in  the 
source  deck.  A listing  of  the  elements  contained  on  these  cards  is 
shown  in  Figure  I.1U. 

Output  Format  - PLM0D2  produces  two  types  of  output  data.  A 7- 
track  magnetic  tape  containing  machine  instructions  for  the  off-line 
plotter  is  generated  during  execution  of  the  program.  The  plots 
ultimately  produced  from  this  tape  have  the  format  of  those  in 
Figure  8,  of  Section  V.  As  an  adjunct  to  these  waveform  plots 
PLM0D2  lists  the  voltage  arrays  that  correspond  to  the  N discrete 
points  plotted,  where  N is  the  number  of  time  increments.  An 
example  of  this  listing  is  shown  in  Figure  I. 11. 

Part  4:  Operational  Stens  - The  chart  shown  in  Figure  1.12 

briefly  describes  the  sequence  of  steps  required  to  exercise  the 
program  as  written  on  a computer  facility  equivalent  to  the 
author's. 
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Figure  1.9  FLOW  CHART  OF  PLM0D2 
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CALL  SYMBOL l *<.50,*,.->,rt.  1«,«  ■*  AHt  OT  AL  Nj  1 t t OF  *.  LL  ^ * *0.7*  -> 

CALL  NUMHf.  R(  999.0*9  99.2,0.  )4*CT  l IS  ,C  . 0 , -1  I 
CALL  SYH30U3.  or,  5.0,0.  ? 1 , 1 JHS  :R r AL  ,0.9,12) 

CALL  SYMBOL ( 999.0 • 999 . 0,  ( .2 ! , SrAl M *0. - , ® ) 

CALL  SYMBOL  < 3.  lb, ->4,5,0.  1*1  H.IMJX— 'OS  1 M Jl. A T l 0 , , " . 0 ♦ 1 M 
X T C t . L L * X M T R 
XTV  'Lt  *RccVL  7 

CALL  SYMBOLIC. C*-5.0,0.W*17Hll'J-  * XMT.ICCLL  , . C , 1 ? ) 

CALL  NUM4EHI99  ,.0,V'-W.0,O.IA,xTC  H.^.O*-!) 

CALL  SYMBOL  ( 999.0 ,999.0,0.  1 4 * 7H  » VOLT  ,0.',7) 

CALL  NOMJ=Rl  999  *0  * 949 . 0 * 0 • l '* , X T V 3LT , .0,-1 ) 

CALL  SYMBOL ( 999. 0,999. 0,7.1  A, 1H) ,0.0, 1 ) 

RCCcLL<CLL  (N) 

RC VOLT  *VLT ( \ ,M) 

CALL  SYM«30U6.3t-5.'»,0.14t17^0A^H  * >CV  < ICUL  ♦''.'tli) 

CALL  NUMBER t 799 .0 * 999 .3 • C • 14 * RCC ELL *G .0* -1 ) 

CALL  SYMOULl 309. 0,999. 0,(.1a, 7M,  VOLT  ,^.0,7) 

CALL  NUM3CR!  999.0  ,9*? .0 * 0 . 14. RCVOLT *n .0  * -l I 
CALL  SYM RJL(99V, 0,099, 0,0.14, 11), C.7 ,1) 

CALL  PLOT ( C . ,-4.,-3) 

CALL  L I NS  t T , XM  TRV  * N T t NC  T ,1 *6,0) 

CALL  OASHL < T .PLOT V , NT  1 NCT , 1 ,0,0) 

CALL  PLOTiO. , 4 * , — 3 1 

if(t:st.le.dcl)  go  to  75 

CALL  SYM8ULr3.4t'». 7*0.?l , 1 *UALTcR£0  Y-AX1  S,  0.  o,  14  ) 

CALL  SYMBOL! 3.6,2.2*0.21 1 1CM3CAL6  FAC rJ  : , r.* , l ’) 

7b  C0NT1\J= 

80  CONTINUE 

CALL  PLOT ( 14.0*0.0* '*9  9 ) 

C PRINT  THE  VOLTAGE  ARRAYS 

VKITr  IS*  1000  )SCRIAL,O0U0 

noc  FJRMATIM  SERIAL  NJ.  *,lA4t*  OAT.  I>  * , 2 A* ) 

WRITMb,  1002  )TM,NT1NCT,NN 

1^02  FORMAT!  *0  PERIOD*1  , El  0 .4 , * 4 Of  T11S  INC  <"MEN  TS  * * , ; C , 

1*  ft  Jr  ODD  HARMONICS  CONS  lUtrKEO*  * *131 
VR 1 T£ (6* 1C00 JREFVlT 

100®  FORMAT! *0  REFER SNCE  VOLTAGE  IN  XMT 1 CELL  IS  # SIC) 

WR1 TE (6, 10091 


LISTING  OP  PLH0D2  (Con't) 
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*JK  I K ^ 

ons 

0H6 

0137 

one 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0143 

0149 

0150 

0151 

0152 

0153 


IV  G LEVEL  31  MAIN  OAT  p * 74277  09/06/57  PAGE  0004 

1009  FORMAT ( / » 2 5X  » * ARRAY  ELEMENTS  ARE  ORDERED  FROM  LEFT  TO  RIGHT*. 

1/j  25X » • AMD  TOP  TO  BOTTOM  OF  THU  TABLES* I 

00  200  M-l.CELLKT.l 
WRITE (6. 1010  1C  ELL IN  I 
101C  FORMAT ( • j CELL  ' , 121 
HQ  190  M-1,9.1 
IC(VLT(M,M1 .EO.OIGO  TO  190 
VOLT* VLT IN. Ml 
WRITE ( o » 101 2 1 VOLT 
1012  FORMAT! *0  VOLTAGE*. I2J 

WRITE  !6,  ICCO) ! V!  I .VLTM.M)  »CrLL<  Ml  1 , I“l  .NTINCT,  1 ) 

1020  FORMAT! * *.l OF  11*71 
190  CONTINUE 
20C  CONTINUE 

c 

C 9RINT  THE  REFERENCE  VOLTAGb  ARRAY 

r 

WRITE  15.  J01O)XMr* 
m*rtfvlt 

WRITE (o.IOIC)M 

WRITE (6. 1020  MVII .M, YMTR) ,1*1 .NTIMCT, 1) 

STOP 

END 


LISTING  OF  PLK0D2  (Con’t) 
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IA  -45 ,079 


DATA 

CARD 

NUMBER 

VARIABLE 

NAME 

DESCRIPTION 

CARD 

COLUMNS 

FORMAT 

COMMENTS 

1 

NTINCT 

NUMBER  OF  TIME 
INCREMENTS  PER 
PULSE  PERIOD 

1 -3 

13 

DETERMINES  GRANULARITY 

of  plots  proouceo 

* 

2 

refvlt 

NUMBER  OF  VOLTAGE 
IN  TRANSMITTER  CELL 
WHICH  IS  CONSIDERED 

1 - 2 

12 

the  REFERENCE  VOLTAGE 
IS  PLOTTED  ON  EVERY 

THE  REFERENCE. 

(•  Q 4 * V (4) , 7 * V{  7 ) , 
ETC.) 

PLOT  TO  GIVE  A POINT 
OF  COMPARISON  FOR 

both  amplitude 

AND  PHASE 

3 

CELL { N ) 

NUMBER  OF  CELL 
CONTAINING  the 
VOLTAGES  TO  BE 

1 - 2 

12 

CELL  NUMBERS  MUST 
BE  RIGHT  JUSTIFIED 

plutted 

VLT  { N,M ) 

THE  NUMBER  OF  THE 
VOLTAGE  IN  CELL  (N) 
TO  BE  PLOTTED 

3 

I ! 

the  numbers  of  the 

VOLTAGES  IN  THIS  CELL 

4 

1 

I 

1 

THAT  ARE  TO  BE 
PLOTTED  ARE  ENTERED 
SEQUENTIALLY  WITHOUT 

IMBEDDED  BLANKS 
OR  DELIMITERS 

T 
1 1 

XI 

/*.g.  I2JLT2.  * Plot  \ 

( VOLTAGES  V{3),V{7),&  V(2)  ] 

\0F  CELL  # 12  J 

(K  £34) 


EACH  CARO  HAS  THE  SAME 

FORMAT  AS  # 3 ABOVE.  A FIXED 
NUMBER  OF  DATA  CARDS  IS  NOT  REQUIRED, 
BUT  ONE  CARD  IS  NEEDED  FOR  EACH 
CELL  CONTAINING  A VOLTAGE  TO  BE 

plotted 


A NUMBER  >50  IS 
PLACED  HERE  ASA 

TRAILER  TO  signify 
THE  END  OF  THE 


I - 2 


12 


CUSTOMARILY  : 99 


DESI  RE  D 
LISTI  NG 


PLOT 


Figure  T.IO  PLM0D2  INPUT  DATA  CARD  FORMAT 
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i mM 


c c 
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FIGURE  I. II 


A -45 ,075 


Figure  1. 12  OPERATIONAL  FLOW  CHART 
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